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Abstract 

An angular analysis of the decay —)• K*^K*^ is performed using pp collisions corre¬ 
sponding to an integrated luminosity of 1.0 fb“^ collected by the LHCb experiment at 
a centre-of-mass energy ^/s = 7 TeV. A combined angular and mass analysis separates 
six helicity amplitudes and allows the measurement of the longitudinal polarisation 
fraction /z, = 0.201 ± 0.057 (stat.) ± 0.040 (syst.) for the ^ A*(892)°A*(892)° 
decay. A large scalar contribution from the Aq( 1430) and Aq( 800) resonances is 
found, allowing the determination of additional CP asymmetries. Triple product and 
direct CP asymmetries are determined to be compatible with the Standard Model 
expectations. The branching fraction H(H° —A* (892)°A* (892)°) is measured to 
be (10.8 ± 2.1 (stat.) ± 1.4 (syst.) ± 0.6 {fd/fs)) x 10“°. 
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1 Introduction 


The —)■ K*^K*^ decay is mediated by a 6 —)■ sdd flavour-changing neutral current 

(FCNC) transition, which in the Standard Model (SM) proceeds throngh loop diagrams 
at leading order. This decay has been discussed in the literature as a possible held for 
precision tests of the SM predictions, when it is considered in association with its U-spin 
symmetric channel 5° —)■ K*^K*^ In the SM, the expected CP violation in the 

former is very small, (9(A^), with approximate cancellation between the mixing and the 
decay CKM phases |^. When a scalar meson background is allowed, in addition to the 
vector-vector meson states, six independent helicities contribnte |^. 

In this paper, a search for non-SM electroweak amplitudes is reported in the decay 

—)■ K^tt~K~ 7r~^, with Ktt mass close to the iF*(892)° mass, throngh the measnrement 
of all CP-violating observables accessible when the flavonr of the bottom-strange meson 
is not identified. These observables include triple products (TPs) and other CP-odd 
qnantities |^, many of which are, as yet, experimentally nnconstrained. Triple prodncts 
are T-odd observables having the generic strnctnre ■ {v 2 x v^) where vi is the spin or 
momentum of a final-state particle. In vector-vector final states of B mesons they take 
the form q ■ (ei x € 2 ) where q is the momentnm of one of the final vector mesons and ei 
and €2 are their respective polarisations. Triple products are also meaningful when one of 
the final particles is a scalar meson. 

Theoretical predictions based on perturbative QCD for the decay of B mesons into 
scalar-vector final states iPQ(1430)P'*(892)° have been recently investigated, yielding 
branching fractions comparable to those of vector-vector final states [^, which have been 
previously available 0- The P° —)■ K*^K*^ decay was first observed with 35 pb ^ of LHCb 
data reporting the measurement of the branching fraction and an angular analysis. 
A remarkably low longitndinal polarisation fraction was observed, compatible with that 
found for the similar decay P° —)■ [^, and at variance with that observed in the mirror 
channel B^ —)■ K*^K*^ 10 and with some predictions from QCD factorisation [7,11. 

An updated analysis of the P° —?• K^7i~K~n~^ final state is reported in this publication, 
in the mass window of ±150MeV/c^ around the iC*(892)° (hereafter referred to as K*^) 
mass for K^tt~ and K~'k^ pairs. A description of the CP observables is provided in Section 
1^ the LHCb apparatus is summarised in Section]^ and the data sample described in Section 
1^ Triple products and direct CP asymmetries are determined in Section A measurement 
of the various amplitudes contributing to B^ —)■ K^tt~K~'k^ is performed in Section]^ 
under the assumption of CP conservation. These include the polarisation fractions for 
the vector-vector mode B^ —)■ In light of these resnlts, the measnrement of the 

branching fraction i3(i?° —)■ is updated in Section]^ Conclusions are summarised 

in SectionThese studies are performed using 1.0 fb“^ of pp collision data from the LHC 
at a centre-of-mass energy of y/s = 7TeV and recorded with the LHCb detector. 
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2 Analysis strategy 


Considering only the S-wave {Ji ,2 = 0) and P-wave (Ji ,2 = 1) prodnction of the Kn 
pairs, with Ji 2 the angnlar momentnm of the respective Ktt combination, the decay 
—)■ 'K~)j^{K~Ti ^)can be described in terms of six helicity decay amplitudes. A 

two-dimensional ht to the K^tt~ and K~tt^ mass spectra, for masses up to the ib}(1430)‘’ 
resonances, hnds a small contribution (< 1%) of tensor amplitudes when projected onto 
the Ktt mass interval used in this analysis, and thus these amplitudes are not considered. 
Three of the above amplitudes describe the decay into two vector mesons, commonly 
referred to as P-wave amplitudes, —)■ V 1 V 2 with Vi = and V 2 = A'*®, with the 

subsequent two-body strong-interaction decay of each of the vector mesons into a Ktt pair. 
Each amplitude corresponds to a different helicity (L^ = 0, -|-1, —1) of the vector mesons 
in the hnal state with respect to their relative momentum direction, Hq, and if_. It is 
useful to write the decay rate in terms of the amplitudes in the transversity basis, 

Ao = i/o, A\\ = ^{H+ + H_) and = (1) 


since, unlike the helicity amplitudes, they correspond to states with dehnite CP eigenvalues 
(^11 =^0 = 1 and r]_[_ = — f). The P-wave amplitudes are assumed to have a relativistic 
Breit-Wigner dependence on the Ktt invariant mass. 

In addition, contributions arising from decays into scalar resonances or non-resonant 
Ktt pairs need to be taken into account within the mass window indicated above. The 
amplitudes describing this S-wave configuration are Ays, Asv and corresponding to 
the following decay^ 

Ay^ : 

A 5 v : 50^ (ib+ 7 r-)oi?*o and ( 2 ) 

A55 : i?°^(ir+vr-)o(ib- 7 r+)o. 


where the subscript denotes the relative orbital angular momentum, J, of the pair. The 
scalar combinations {Ktt)q are described by a superposition of a broad low-mass structure 
related to the A’q( 800) resonance 12 and a component describing the A’q( 1430) resonance. 

Unlike the K*^K*^ hnal state, the S-wave conhgurations SV and VS dehned in 
Eq. ([^ do not correspond to CP eigenstates. However, one may consider the following 
superpositions 


|s+) = ^(\K’\K-7r*)„) + \(K*7r-)oK’'>)) and 

la-) = ^(|A'-»(ii-V)„)-|(ii'+7r-)„F*»)), (3) 

which are indeed CP eigenstates with opposite CP parities {7]s+ = — 1 and rjg- = -|-1). 
Therefore, it is possible to write the full decay amplitude in terms of CP-odd and CP-even 

^Note that both 5° and can decay into these final states. 
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( 4 ) 


amplitudes (the SS final configuration is a CP eigenstate with rjss = 1) by defining 


i^vs + ^sv) and 


= {^vs - Asv) 


2.1 Angular distribution 

The angles describing the decay, hi = {9i,92,^p}, are shown in Fig. where ^ 1 ( 2 ) is the 
angle between the direction of K~^C) meson and the direction opposite to the i?-meson 
momentum in the rest frame of Fi( 2 ) and (p is the angle between the decay planes of the 
two vector mesons in the i?-meson rest frame. In this angular basis, the differential decay 


rate describing this process is expressed as 13 

d^r 


dhl dmi dm 2 dt 


= N 


Ao(t) cos 6*1 cos 6*2 + 


A||(t) 

72 


sin 6*1 sin 02 cos ip 


+i ^ sin 9i sin 02 sin p ) Ali(mi)Ali(m 2 ) 

V 2 

^tjt) ( 

76 V 


76 


cos0iAli(mi)Alo(’^2) — cos92Mo{mi)M.i{m2) 
cos9iMi{rrii)Mo{m2) + cos92Mo{rrii)Mi{m2) 


■'^Mo{mi)Mo{m2) 


(5) 


where the different dependences of P-wave and S-wave amplitudes on the two-body masses 
mi = M{K^tt~) and m 2 = have been made explicit in terms of the mass 

propagators and N is an overall normalisation constant. The time evolution 

induced by mixing is encoded in the time-dependence of the amplitudes Ak{t) 

{k = 0, II, T, s+, s“, ss) 

Ak{t) = g+{t)Ak + g-{t)Ak, (6) 

where Ak = Akit = 0 ) and the time-dependent functions g±{t) are given by 





(7) 


with Vi H {jn^^u) being the width (mass) of the B^ light (L) and heavy (FT) mass 
eigenstates. 

The decay rate of the CP-conjugated process, P° —)■ (iF“ 7 r+)(iF+ 7 r“), can be obtained 
by exchanging each amplitude Ak by rjkAk, where gk is the CP eigenvalue of the final state 
described by A^ Q. In this paper, due to the limited size of the available data sample, 
no attempt is made to identify the flavour of the initial meson at production, thus 
suppressing the sensitivity to direct and mixing-induced CP asymmetries. Nevertheless, 
CP violation can still be studied through the measurement of triple product asymmetries 
and S-wave-induced direct CP asymmetries. 
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Figure 1: Definition of the angles involved in the analysis of —)> K*^K*^ decays. 


2.2 Triple product asymmetries 


Two TPs can be defined in B meson decays into pairs of vector particles |^[5, 14 

Ti = (riyj X hv^) ■ Pvi = sin93 and 


( 8 ) 


T 2 = 2{hvi ■ nv 2 ){nv^ X nyj • pv^ = sin2(p, 


(9) 


where hy. (i = 1, 2) is a nnit vector perpendicular to the Vi decay plane and pv^ is a unit 
vector in the direction of Vi in the rest frame. The observable asymmetries associated 


with these TPs can be calculated from integrations of the differential decay rate as 14 


a^(t) 


a^(t) 


P(cos6'i cos 02 sin (p > 0,t) — P(cos0i cos 02 sin (p < 0,t) ^ 

P(cos0i cos 02 sin > 0,t) + P(cos0i cos 02 sin < 0,t) 

P(sin2(p > 0,t) — P(sin2(p < 0,t) 

P(sin2(p > 0,t) + P(sin2(p < 0,t) 


( 10 ) 

( 11 ) 


Nonzero TP asymmetries appear either due to a T-violating phase or a T-conserving phase 
in conjunction with hnal-state interactions. 

When these asymmetries are measured in a sample where the production flavour is not 
identihed, they become “true” CP-violating asymmetries, assuming that CPT is conserved. 
This is manifest when and are written in terms of the amplitudes dehning the decay 
rate in Eq. (|^, 

2 p2u 1 

= - —Xm{A 2 _A*Q) and (12) 

TT u 

al = --^XmiAM, (13) 

TT U " 


with T> = \Ao\'^ + |.A||p -h \A±\‘^ + -h \A~\‘^ + Taking into account these 

expressions and that the decay rate contribution associated with the CP-odd amplitude 
A±_ changes sign under the CP transformation, the asymmetries measured in the untagged 
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sample, are proportional to the CP-violating interference terms Xm{A^AQ y — ^±v4,q | 
Using Eq. these terms can be written as 


Xm(^_L^Qyi — Xl_L^o,||) = 2^ 


— Cat 




cosh 

sinh 


Ar, 

2 

AE. 


(14) 


where AE^ = E^ — E^^, E^ = (E^, + Ej:^)/2 and 0mix is the phase in P°-P° mixing. The 
coefficients Xm^A^A^ y — A^A^ y) and Xm[[A^AQ y -|- are TP and mixing- 

induced TP asymmetries, respectively, and are CP-violating quantities j^. In the analysis 
presented in this paper, only the time-integrated asymmetries 

Alp = —f Xm{A±A*Q) dt and 
TT D J 

Ap = —— [ Xm{A±Add dt, 

n U J 

are measured {D = J Vdt), with no identification of initial P° flavour. Thus CP-violating 
linear combinations of the above observables are accessible. 

When the S-wave contribution is taken into account, two additional CP-even amplitudes, 
A~ and Xlss, interfere with XI_l, and give rise to two additional CP-violating terms. Further 
asymmetric integrations of the decay rate, analogous to those in , lead to the following 
observables 


(15) 

(16) 



E((cos 6 *i + cos 6 * 2 ) sin(^ > 0) — E((cos 6 *i -|- cos 6 * 2 ) sin(^ < 0) 
E((cos 61 + cos 62 ) sin 99 > 0) -|- E((cos 6*1 + cos 6 * 2 ) sin (^ < 0) 




A±As *) Mi{m)Ml{m)) dm 


(17) 


and 


^4 _ r(sin(^ > 0 ) — r(sin(^ < 0 ) 

E(sin(^ > 0) + E(sin(^ < 0) 


dvr 1 


Xm Mi{m)M*Q{rn)) dm. 


(18) 


where the mass integration extends over the chosen Kn mass window. It is performed 
over the product of mass propagators of different resonances, times specific CP-vioIating 
observables involving Xlj,. 

Since is also CP-odd, its interference terms with the CP-even amplitudes change 
sign under P° to interchange. Consequently, four new CP-violating asymmetries are 
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accessible from —>■ i^+vr K decays, 

_ r(cos 6 'i cos 6 ' 2 (cos 6 'i — COS 62 ) > 0) — r(cos 6 'i cos 6 ' 2 (cos— cos 6 ' 2 ) < 0) 
^ r(cOs 6 'i cos 6 ' 2 (cOS 01 — cos 02) > 0) + r(cos0i COS 02 (cOS 0 i — cos 02 ) < 0) 

/H 1 r r 


1 

5V^V 


9 hze 


~ AQ)Mo{m)M\{rn)) dm 
+5 j ne{{AtAl,-AtAl,)Mi{m)Ml{m))<lm 


42 


r((cos0i — cos 02 ) COS</9 > 


, 0) — r((cOS0i — COS 02) COS</9 < 0) 

r((cOS0i — COS 02) COS</9 > 0) + r((cOS0i — COS 02) COS < 0) 
QO 1 /“ 

y - ^Mi)-Mo(m)>lKH) dm, 


43 


r((cos 01 

r((cos0i 

2 V 2 1 r 

5^3 L 


+5 


cos 02 ) > 0 ) — r((cos 0 i — cos 02 ) < 0 ) 
cos 02 ) > 0 ) + r((cos 0 i — cos 02 ) < 0 ) 

TZe ((^j '^0 ~ A'^Al)Mo{m)Ml{m)) dm 
[ -^+^*JAfi(m)Afo(m))dm 


( 19 ) 


( 20 ) 


( 21 ) 


and 




r((cOS^ 01 — COS^ 02) > 0) — r((cOS^ 01 — COS^ 02) < 0) 

r((cOs 2 01 — COS^ 62 ) > 0 ) + r((cOs 2 01 — COS^ 02) < 0 ) 


= pKe {AtA, 


— * 
s 


-4J-4;*) ■ 


( 22 ) 


Some of these terms have the form 7le{A'^Al — A'^AX), with A; = 0, ||,s“,ss, which is 
characteristic of direct CP asymmetries. 

As shown above, TP and several direct CP asymmetries are accessible from untagged 
Bg —)■ K*^K*^ decays, provided that a scalar Kti background component is present. These 
CP-violating observables are sensitive to the contributions of FCNC processes induced 
by neutral scalars, which are present, for example, in models with an extended Higgs 
sector. Constraints on possible FCNC couplings of Higgs scalars have been recently 


examined 15 T 6 


Non-zero values of or would allow the characterisation of those 

operators contributing to the effective Hamiltonian. In particular an enhanced contribution 
from A^^’^ with respect to the other observables, would reveal stronger [V — A) x {V + A) 
{LR) and (H -|- A) x (H — A) {RL) components with respect to RR and LL operators in 
the above models [^. 
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2.3 Angular analysis 

Assuming that no CP violation arises in this decay, an angular analysis of the decay products 
determines the polarisation fractions of the 5° —?■ K*^K*^ decay and the contribution of 
the various S-wave amplitudes. The time-integrated decay rate can be expressed as 


d^r 

did dmi dm2 


21 

N E A'„(mi,m2)F„(ld), 

n=l 


(23) 


where the functions Kn contain the dependence on the amplitndes entering the decay, 
with their corresponding mass propagators, and N is an overall normalisation constant. 
The Kn functions are given in Table together with the decay angle fnnctions All 
terms proportional to the TP and S-wave-induced CP asymmetries (n = 5, 6, 21 and the 
symmetric A+ <H- A~ terms in n = 8 — 11,13 — 16) cancel nnder the assnmption of CP 
conservation. 

The dependence of each amplitude on the invariant mass of the K~^7r~ and iP“7r+ pairs 
is given by the propagators A4.j{m) oc TZj{m) x m/q, where q is the momentnm of each 
meson in the rest frame of the Kir pair 


q = 


VK - {Mn + MKf){m^ - {M^ - MkT) 


2 m 


(24) 


The P-wave propagator, J = 1, is parameterised using a spin-1 relativistic Breit-Wigner 
resonance function 

^ / X MiPi(m) , , 

^ (M2-m2)-zMiPi(m)' 

The mass-dependent width is given by 

3 


Pi(m) = P 


Ml 1 -|- Cq( f q 


m 1 -|- Cq^ \qi 


(26) 


where Mi and Pi are the iP*'^(892) resonance mass and width, r is the interaction radins 
and qi corresponds to Eq. (24) evaluated at the resonance position (Mi). 


To describe the S-wave propagator, the LASS parameterisation is used, 

which is an effective-range elastic scattering amplitude, interfering with the iPg(1430) 
resonance. 


7lo{m) oc 


1 


- + e 


2 iS/^ 


Moro(m) 


where 


cot Si 3 — i 

Po(m) = To— ( - 


M2 - m2 -iMoPo(m)’ 
Mo f q 


m \qo 

and the non-resonant component is described as 

cot d /3 =-h -bq. 

aq 2 


(27) 

(28) 

( 29 ) 
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Table 1: Untagged time-integrated terms used in the analysis, under the assumption of no CP 
violation. 


n 

Kn 

Fn 

1 

^\Ao\^\Mi{mi)\'^\Mi{rn2)\'^ 

COS^ 01 COS^ 02 

2 

I^|A||nMi(mi)p|Mi(m2)P 

\ sin^ 01 sin^ 62 cos^ ip 

3 

r^|A^nMi(mi)nMi(m2)P 

1 sin^ 01 sin^ 02 sin^ ip 

4 

j^^|A||||ylo| cos 5|||Ali(mi)p|A4i(m2)P 

-Aj= sin 201 sin 202 cos ip 

5 

0 

— -Aj= sin 201 sin 202 sin ip 

6 

0 

— \ sin^ 01 sin^ 02 sin 2ip 

7 

2 ('rJ +'rj )\Mi{mi)\^\Mo{m2)\‘^ 

I cos^ 01 

8 

A4^(m2)2Vlo(m2))|Ali(mi)p 

— ^ cos^ 01 cos 02 

9 

-^^\A~\\A\\\ne{P^^‘’ M\{m2)MQ{m2))\Mi{rni)\^ 

— sin 201 sin 02 cos ip 

10 

^ 1 y4+11 yl_L ) AIq (m 2 ) Alo (m- 2 )) 1 All (mi ) p 

sin 201 sin O2 sin ip 

11 

-^^\A-\\Ass\ne{P^^‘> Ml{mi)Mi{mi))\MQ{m2)\‘^ 


12 

2 ('rJ +'rj )\Aio{mi)\^\Mi{m2)\‘^ 

1 COS^ 02 

13 

-^^\A~\\AQ\ne{P^‘‘ Al^(mi)Alo(w-i))|-Adi(m2)P 

COS 01 COS^ 02 

14 

- ^ lY 1^711^11 -'^ 11 ) All (mi) Alo(mi)) 1 All (m 2 ) P 

^ sin 01 sin 202 cos ip 

15 

^j^|A+||A_L|Xm(e*('^-L-'5^)Alo(mi)Alo(mi))|Ali(m2)P 

— ^ sin 01 sin 202 sin ip 

16 

-^^\A-\\Ass\Re{P^^‘> Ml{m2)Mi{m2))\Mo{rni)\'^ 


17 

('p"J ' p'J )7^e(Al^(ml)Alo(?7^2)7Vlo(ml)All(m2)) 

— 1 cos 01 cos 02 


18 ’^\Ass\^\Mo{mi)\^\Mo{m2)\‘^ g 

19 ^\Ass\\AQ\Pe{P^’^‘‘M\{mi)M\{m2)Mo{mi)Mo{rn2)) -|cos0icos02 

20 M\{mi)M\{m 2 )M.o{rni)M.o{m 2 )) — ^ sin^i sin 02 cosLp 


^ sin 01 sin 62 sin ip 


21 
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Table 2 : Parameters of the mass propagators used in the fit. 



{Ki 
J = 0 

A' 

17 

18 


AT* (89^ 

j = i [5| 

Mj (MeV/c2) 

1435 ± 5 d 

Z 5 

895.81 ±0.19 

Tj (MeV/c^) 

279 ± 6 ± 21 

47.4 ±0.6 

r (GeV-i) 

- 




3.0 ±0.5 

a ( GeV^) 

1.95 ±0.09 ±0.06 

- 

b (GeV^) 

1.76 ±0.36 ±0.67 

- 


The values of the mass propagator parameters, including the resonance masses and widths, 
Mj and Tj, and the the scattering length (a) and effective range (6), are summarized in 
TableOther shapes modelling the S-wave propagator, including an explicit Breit-Wigner 
contribution for the iTQ(800) resonance, are considered in the systematic uncertainties. 

The normalisation of the mass propagators 


lAdopdm = / |A^ipdm = l 


(30) 


in the mass range considered, together with the normalisation condition 

|/l„p + |/l||p + + \A;\^ + = 1, 


(31) 


guarantees the definition of the squared amplitudes as fractions of different partial waves. 
The polarisation fractions for the vector mode, —)• are defined as 


/l,II,± — 


l"4o,||,J 


|Ao|2 + |A|||2 + |yl^|2- 

The overall phase of the propagators is defined such that 

arg[A^o(dTi)] = arg[A^i(Mi)] = 0 


(32) 


(33) 


and the convention 5o = arg(74o) = 0 is adopted. Therefore (5||, 5^, 6 ~, 6^ and 6 ss are 
defined as the phase difference between the corresponding amplitude and Aq at the K*^ 
mass pole. As a consequence of the lack of initial 5° or flavour information, the phases 
5 _l and can not be measured independently, and only their difference is accessible to 
this analysis. 


3 The LHCb detector 

The LHCb detector |20pl] is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < t] < 5, designed for the study of particles containing 6 or c quarks. The 
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detector includes a high-precision tracking system consisting of a silicon-strip vertex 
detector surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4 Tm, and three stations 
of silicon-strip detectors and straw drift tubes placed downstream of the magnet. The 
tracking system provides a measurement of momentum, p, of charged particles with a 
relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200GeV/c. The 
minimum distance of a track to a primary vertex, the impact parameter (IP), is measured 
with a resolution of (15 -|- 29/pt) M-ni, where px is the component of the momentum 
transverse to the beam, in GeV/c. Different types of charged hadrons are distinguished 
using information from two ring-imaging Cherenkov detectors. Photons, electrons and 
hadrons are identihed by a calorimeter system consisting of scintillating-pad and preshower 
detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons are identihed 
by a system composed of alternating layers of iron and multiwire proportional chambers. 
The online event selection is performed by a trigger, which consists of a hardware stage, 
based on information from the calorimeter and muon systems, followed by a software stage, 
which applies a full event reconstruction. 

In the analysis presented here, all hardware triggers are used. The software trigger 
requires a multi-track secondary vertex with a signihcant displacement from the primary 
pp interaction vertices (PVs). At least one charged particle must have a transverse 
momentum px > 1.7 GeV/c and be inconsistent with originating from a PV. A multivariate 
algorithm 


22 identihes secondary vertices consistent with the decay of a 6 hadron. 


Simulated —>■ K*^K*^ events are used to characterise the detector response to signal 

events. In the simulation, pp collisions are generated using Pythia 23 with a specihc 


LHGb conhguration 24 . Decays of hadronic particles are described by EvtGen [25] , 
in which hnal-state radiation is generated using Photos The interaction of the 


generated particles with the detector and its response are implemented using the Geant 4 
toolkit 


27 as described in Ref. 28 


4 Event selection and signal yield 


The event selection is similar to that used in the previous analysis j^. K*^ candidates 
are formed from two high-quality oppositely charged tracks identified as a kaon and pion, 
respectively. They are selected to have px > 500MeV/c and to be displaced from any PV. 
The K~^7t~ and K~n^ pairs are required to have invariant mass within ±150MeV/c^ of the 
known K*^ mass, which corresponds to 74% of the total phase-space for —)■ 

and Px > 900MeV/c. Each B^ candidate is constructed by combining a K*^ and K*^, 
requiring the four tracks to form a good vertex well-separated from any PV. The B^ 
candidate invariant mass is restricted to be within the interval [5100,5866] MeV/c^ and its 
momentum vector is required to point towards one PV. 

In order to further discriminate the B^ —)■ K*^K*^ signal from the combinatorial 
background, different properties of the decay are combined into a multivariate discrimina¬ 


tor 29 . The variables combined in the discriminator are the 5° candidate IP with respect 
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to the associated PV, its lifetime and px, the minimum Xip of the four daughter tracks 
(defined as the difference between the of a PV formed with and without the particle in 
question) with respect to the same PV and the distance of closest approach between the 
two K*^ candidates. The discriminator is trained using simulated —)• K*^K*^ events 
for signal and a small data sample excluded from the rest of the analysis as background. 
The optimal discriminator requirement is determined by maximising the hgure of merit 
Ns /a/ Ns + Nb in a test sample containing signal (S') and background {B) events of the 
same nature as those used in the training sample. 

Differences in the log-likelihood for various particle identification hypotheses (A In Ca-b) 
are used to minimise the contamination from specihc B decays. Contributions from 
B^ —)■ pK*^ and 5° —)■ (f)K*^ modes are reduced by the A In Ck-tt requirements of kaons 
and pions. A small contamination from Aj) —)■ pTi~K~n^ decays is observed and suppressed 
with A In Cp-K requirements. 

An extended unbinned maximum likelihood ht to the mass spectrum of the selected 
B^g —)■ K^tt~K~'k^ candidates is performed. The signal is modelled by a sum of two 
Crystal Ball distributions that share common mean and width. The same distribution 
is used to describe the B^ decay into the same hnal state. Components for B^ —)■ (pK*^ 
and A° —)■ pn~K~n^ decays are included in the fit with shapes extracted from simulated 
events. The contribution from B^ —)■ pK*^ decays is estimated to be negligible from 
simulation studies. Finally, partially reconstructed B decays are parameterised using 
an ARGUS distribution and the remaining combinatorial background is modelled 
using an exponential function. The ht result is shown in Fig. A total of 697 ± 31 
B^ —)■ K^7t~K~tt^ decays is obtained. 

4.1 Acceptance properties 

Effects introduced in data due to the geometry of the detector and to the selection 
requirements need to be taken into account in the measurement. 

The study of simulated B^ —)■ K*^K*^ events shows that the detection and selection 
efficiency is not uniform as a function of the decay angles 9i and 6 ^ 2 , but has no dependence, 
at the level of precision needed for this analysis, on p and on the invariant mass of the 
two Ktt pairs, rrti and m 2 - The acceptance decreases as cos^j approaches 1. This feature 
is mainly induced by the requirement on the minimum px of the daughter pions. This 
effect is modelled by a two-dimensional function in cos 6 *i and cos 6 * 2 , which is extracted 
from simulation. 

Since the trigger system uses the px of the charged particles, the acceptance effect 
is different for events where signal tracks were involved in the trigger decision (called 
trigger-on-signal or TOS throughout) and those where the trigger decision was made using 
information from the rest of the event (non-TOS). The data set is split according to these 
two categories and a different acceptance correction is applied to each subset. 


11 




5 


_5^-1-,-,-,-1-,-,-,-1-,-,-,-1-J 

5200 5400 5600 5800 

M(KVK“;i^) [MeV/c^] 


Figure 2: Invariant mass distribution for selected K'^'n~ candidates. The (blue) solid 
line is the result of the fit explained in the text. The and signal peaks are shown as 
as dashed-dotted lines (pink and dark green, respectively). The various peaking background 
components are represented as dotted lines: (red) B^ —)• (pK*^, (green) —?■ p7r~K~ tt^ and 
(light blue) partially reconstructed decays. The (grey) dotted line is the combinatorial background 
component. The normalised residual (pull) is shown below. 


5 Triple product and direct CP asymmetries 


Triple products and direct CP asymmetries are calculated for —)■ if'+vr K 7r+ using 
Eqs. (10) and after time integration, and Eqs. (|Tg-(@ from those candidates with 
a four-body invariant mass within ±30 MeV/c^ of the known mass. The background 
in this interval, which is purely combinatorial, is subtracted according to the fraction 
calculated from the result of the invariant mass ht, /bkg = (3.44 ± 0.34)%. The angular 
distributions of the background are extracted from the upper mass sideband, dehned by 
MTT~K~Ti^) > 5550 MeV/c^. Acceptance effects are then corrected in the signal 
angular distributions. The measured asymmetries are listed in TableFrom the definitions 
given in Sect. 2.2, correlations of the order of 5% are expected among these asymmetries. 


with the exception of A\, and where the correlation is calculated to be close to 90%. 

The main systematic uncertainty in these measurements is associated to the angular 
acceptance correction. Discrepancies in the pt spectra and the particle identihcation 
efficiencies between data and simulation are used to modify the acceptance function 
obtained from simulation. Systematic uncertainties are determined from the variation in 
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Table 3: Triple product and direct CP asymmetries measured in this analysis. The first 
uncertainties are statistical and the second systematic. 


Asymmetry 

Value 




0.003 

± 

0.041 

± 

0.009 

Al 

0.009 

± 

0.041 

± 

0.009 

A^ 

0.019 

± 

0.041 

± 

0.008 

A^ 

-0.040 

± 

0.041 

± 

0.008 


-0.061 

± 

0.041 

± 

0.012 


0.081 

± 

0.041 

± 

0.008 


-0.079 

± 

0.041 

± 

0.023 

Ah 

-0.081 

± 

0.041 

± 

0.010 


the measured asymmetries when this modified acceptance is used. Systematic effects are 
found to be larger in case of the four direct CP asymmetries, in particular for which 
has a strong dependence on cosdi 2. In addition, the lifetime-biasing selection criteria 
have a slightly different effect on the various amplitudes, which correspond to decays with 
different effective lifetimes, due to the width difference between mass eigenstates 
This could induce a bias in the measured TP and direct CP asymmetries. A set of 
simulated experiments is performed to estimate the impact of the lifetime acceptance in 
the eight quantities. The observed deviations are small and are used to assign a systematic 
uncertainty. Finally, the effect of the uncertainty in the background contribution is 
estimated by changing the background fraction and the parameters in the background 
model within their statistical uncertainty and recalculating the asymmetries. 
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6 Angular analysis 

The magnitudes and phases of the various amplitudes contributing to the —?■ 
K^tt~K~tt^ decay are determined using a five-dimensional fit to the three helicity angles 
(12) and to the invariant mass of the two Ktt pairs (mi, m2) of all candidates with a 
four-body invariant mass \M{K^, K ~— m^o| < 30MeV/c^. 

The model used to describe the distribution in these five variables is given by 


P{Vt,mi,m2) = (1 - /bkg)i^(f^,"ii,"i2) X £(12) -F /bkg^bkg(f^,"ii, m2), (34) 


where T(12,mi,m2) is the probability density function in Eq. ([^, £(12) is the acceptance 
function modelling the effects introduced by reconstruction, selection and trigger reported 
in Sect. 4.1, and Fbkg(12, mi, m2) describes the distribution of the background extracted 
from the upper mass sideband. The background fraction, /bkg, is obtained from the result 
of the fit to the invariant mass of the candidates. 

Using this model, an unbinned maximun likelihood fit is performed simultaneously for 
TOS and non-TOS —?• K^n~K~7i~^ candidates, where only the acceptance function and 
the background fraction are different between the two samples. The results of the fit are 
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Table 4: Results of the simultaneous fit to —)• K^'k~K~'k^ TOS and non-TOS candidates 
with \M{K~^,TT~, K~,TT^) — m^o\ < 30MeV/c^ (phases are measured in radians). The first 
uncertainties are statistical and the second systematic. 


Parameter 

Valne 

h 

0.201 ±0.057 ±0.040 

f\\ 

0.215 ±0.046 ±0.015 

\At? 

0.114 ±0.037 ±0.023 

\A-? 

0.485 ±0.051 ±0.019 

IA..P 

0.066 ± 0.022 ± 0.007 

<^11 

5.31 ±0.24 ±0.14 

1- 

1 : 

+ 

1.95 ±0.21 ±0.04 

^7 

1.79 ±0.19 ±0.19 

Sss 

1.06 ±0.27 ±0.23 


summarised in Table Figures and show the angular and Kn mass projections of the 
multi-dimensional distributions. To quantitatively demonstrate the interference between 
the different partial waves, a forward-backward asymmetry is defined for K*^ meson as 
Afb = (^F — ^b)/(^f + ^b), where Np {Ns) is the number of mesons emitted with 
positive (negative) cos 6^1, and analogously for the K*^ meson. Their evolution with the Kn 
invariant mass is shown in Fig. as an additonal projection of the ht result. According to 
Eq. (|^ these asymmetries are proportional to the interference term between A~ and Aq. 

Figureshows the likelihood for the longitudinal polarisation fraction fp, where all the 
other parameters are minimised at each point of the curve, depicting parabolic behaviour 
around the minimnm. Additionally, confidence regions in the plane are shown. 

The most important systematic nncertainties in the measnrement of the different 
amplitudes, phases and polarisation fractions are summarised in Table They arise mainly 
from uncertainties in the modelling of the Kn mass distributions and from the assumption 
that the five-dimensional acceptance factorises into a product of two-dimensional functions. 
To better exploit the statistical power of the simulated sample in the less populated regions 
of the phase space, e.g. the tails of the mass distribution, the angular and mass acceptances 
are assumed to factorise. An alternative model is tested that allows for correlation between 
the angnlar distribntion and the Kir invariant mass, nsing a two-dimensional function in 
{cos 6i,mi), nniversal for K*^ and K*^ decays. The ht is repeated with this acceptance 
model and a systematic nncertainty, Uacc , is determined from the variation with respect 
to the nominal ht resnlt. An additional nncertainty acconnts for the limited size of the 
simulated samples, Usim. 

To test the accuracy of the simulation, kinematic distribntions, snch as those of the 
Pt of hnal-state particles, are compared between data and simnlation. Since the inpnt 
amplitudes used in the generators are diherent from those measured in data, an iterative 
method is dehned to disentangle the discrepancies associated with a diherent physical 
distribntion. This procedure supports the quality of the simulation, and allows for the 
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Figure 3: Results of the simultaneous fit to —)• K^7:~K~tt^ candidates (blue solid line) in the 
three helicity angles. The dots represent the data after background subtraction and acceptance 
correction. The red dashed line is the P-wave component, the green dashed line is the S-wave 
component and the light-blue dashed line represents the Aq interference term. 






Figure 4: Projections of the model fitted to —)• K~7r~^ candidates (blue solid line) in 

the (top) invariant mass of RTtf pairs and (bottom) cos 9 asymmetries as functions of Ktt mass 
(mi, m 2 ). The dots represent the data after background subtraction and acceptance correction. 
The red dashed line is the P-wave component, the green dashed line is the S-wave component 
and the light-blue dashed line represents the Aq interference term. 
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Figure 5: (Left) Profile likelihood for the parameter f^. (Right) Regions corresponding to 
Aln£ = 0.5, 2 and 4.5 (39%, 87% and 99% conhdence level) in the \A~\‘^ - fi plane. 


Table 5: Systematic uncertainties in the measurement of the magnitude and phase of the different 
amplitudes contributing to the —)• K^'k~ K~tt^ decay. 


Parameter 

^acc 

^sim 

(Tyw 

^mass 

^res 

Total 

h 

0.031 

0.010 

0.010 

0.021 

0.006 

0.040 

f\\ 

0.008 

0.008 

0.004 

0.005 

0.007 

0.015 

\At? 

0.019 

0.005 

0.002 

0.011 

0.003 

0.023 


0.007 

0.007 

0.010 

0.003 

0.012 

0.019 


0.003 

0.001 

0.000 

0.005 

0.003 

0.007 

(5|| 

0.130 

0.037 

0.042 

0.005 

0.025 

0.144 

ix - 

0.016 

0.019 

0.000 

0.017 

0.027 

0.040 


0.160 

0.036 

0.075 

0.033 

0.030 

0.186 


0.096 

0.076 

0.188 

0.018 

0.044 

0.229 
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determination of the associated systematic uncertainty, Urw 

Several alternative models for the parameterisation of invariant mass propagators are 
used and a systematic uncertainty, cXmass) for the ht parameters is estimated from the 
variation of the ht resnlts. The main contribntion to this nncertainty comes from the 
S-wave mass propagator, which is modelled by the LASS parameterisation in the 
nominal ht. A combination of two spin-0 relativistic Breit-Wigner distributions with the 


mean and width of the ^”0(800) and ^0(1430), respectively 19 , and a single contribution 
from A'q( 1430) are also used. 

Additional small uncertainties are considered to account for the effect of the invariant 
mass resolution, the lifetime acceptance and possible biases induced by the htting method 
(cTres) • 


7 Measurement of 


The branching fraction of the vector mode —)■ K*^K*^ is npdated with respect to the 
previons resnlt |^. This measurement is normalised nsing the —)• (j)K*^ decay, with 

0 —)■ K^K~ and K*^ —)■ K~^7r~, which has a topology similar to the signal decay and a 
well-known branching fraction. 

The selection of B^ —)■ (pK*^ decays is performed snch that it closely resembles the 
selection of B^ —)■ K*^K*^ decays, except for particle identihcation criteria. In particular, 
the requirements related to the vertex dehnition and the kinematic properties of the 
charged particles are identical. Figure shows the invariant mass of the hnal-state particles 
for the selected candidates. 

The ratio of branching fractions for signal and reference decay channels is given by 


B{B^ 


^ ^BO-^4>K*o ^ Xf^{B^ —)■ (pK*^) 

_ trig _ \ , tDO V 

^ A'so X fBO^K*OK*o fd 13{(p ^ K^K~) 

"" Nbo X fBO^^K*o ''Ts'' B{K*^^K+7r-y 


where fd/fs is the ratio of probabilities for a b quark to form a or a B^ meson 34 
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K+ 


The qnantities Nbq and Nbo represent the nnmber of observed candidates for B^ —)■ 
and B^ —)• K^K~ decays, respectively, and are determined from the 


vr 


-A-7r+ 


corresponding fits to the four-body invariant mass spectra. The value of Nbo is reported 
in Sect. The yield Nbo is extracted from an extended unbinned maximum likelihood fit 
to the spectrum in Fig. The B^ signal is modelled by a combination of Crystal Ball 
and Gaussian distributions that share a common mean. Their relative width, fraction and 
parameters describing the tail of the Crystal Ball function are set to the values determined 
from simnlation. The signal from the recently observed decay 5° —)■ (pK* 


^*0 
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is also 

described using this parameterisation. The mass difference between B^ and mesons is 


fixed to the world average value 19 . The partially reconstructed background is modelled 


using an ARGUS distribution with parameters free to vary in the fit. The combinatorial 
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Figure 6: Invariant mass of the selected K'^K ~combinations and result of the fit to the 
data. The points represent the data and the (blue) solid line is the ht model. The Bg and 
signal peaks are shown as as dashed-dotted lines (pink and dark green, respectively). The 
contribution from partially reconstructed decays is represented as a (light blue) dashed line. The 
(grey) dotted line is the combinatorial background component. The normalised residual (pull) is 
shown below. 


background is parameterised with a decreasing exponential function. A total of 1049 ± 33 
signal decays for the —)• K ~decay are observed. 

The yield of candidates corresponding to the resonant decays, 5° —)■ K*^K*^ and 
—)■ (f)K*^, is given by the purity factors fBO^K*OK*o The ratio of 

combined reconstruction and selection efficiencies, is calculated using —)■ 

and B^ —)■ (j)K*^ simulated events and validated using data. The inefficiency induced 
by the particle identification requirements is then determined separately using large 
calibration samples. The ratio of trigger efficiencies, is compnted throngh a data- 
driven method 37 . Moreover, the overall efficiency for each channel depends on the 


helicity angle distribution of the final state particles, and is encoded into the factors 
Both the purity and Afactors for the B^ —)■ K*^K*^ decay are calculated from the results 
of the angular analysis. Those corresponding to B^ —>■ (f)K*^ decays are calcnlated from 
Ref. 


38 


With the factors summarised in Table the ratio of branching fractions is determined 
to be 


B{B^ 


= 1.11 ± 0.22(stat.) ± 0.12(syst.) ± 0.06(/rf//,). 


(36) 
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Table 6: Summary of relevant quantities in the B{Bg —)• calculation. The factor 

—)• is defined as Xf^{Bg —>■ K*^K*^)/f^o^x*ox*o, and equivalently for the 

B^ —)■ (pK*^ decay. The first uncertainty is statistical, the second systematic. 


N^o 

Nbo 

^ Bg-»K*og*o 


697 ±31 ±11 
1049 ± 33 ± 7 
0.453 ± 0.059 ± 0.040 
1.30 ±0.17 ±0.07 


Using the average —>■ (f)K*^) = (9.73 ±0.72) x 10 ® from the BaBar 18 

measurement^ corrected to take into account different rates of B^B 


and Belle 39 
and B^B^ pair 

production from T(4S') using T{B^B~)/T{B^B^) = 1.055 ±0.025 |^, the result obtained 

is 

^ = (10.8 ± 2.1 (stat.) ± 1.4 (syst.) ± 0.6 {fjfs)) x 10"®. 

The main systematic uncertainties considered are related to the invariant mass ht 
used to determine the signal and reference event yields, the angular correction, and the 
determination of the trigger efficiency. To determine the systematic uncertainty associated 
with the number of candidates, the ht is repeated using different models for the signal and 
background components. The largest variation is assigned as a 1.7% systematic uncertainty. 
A 5% uncertainty is attributed to the trigger efficiency, after calibration of the data-driven 
method applied to both channels using fully simulated events. The systematic uncertainty 
associated with the angular correction Ais the result of the propagation of the systematic 
uncertainties evaluated for the parameters measured in the angular analysis (9%). 

This result supersedes the previous measurement , which used a less sophisticated 
estimate of the S-wave contribution. If rescaled to the same S-wave fraction, both results 
are compatible. 

As a result of B^-B^ mixing, the time-integrated havour-averaged branching fraction 
(B) reported here cannot be directly compared with theoretical predictions formulated 
in terms of the decay amplitudes at t = 0 (Bq). The relation between these branching 
fractions is given by 
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B = /af Bo, with /af = 1 - 


Ar, 

2r. 


(/l ± f\\ ± f±) 


(37) 


Using the decay widths measured in Ref. 33 and the polarisation fractions reported here. 


the correction factor is calculated to be /af = 1-015 ± 0.010. 


8 Conclusions 


The decay B^ —)■ iT+vr K is studied using pp collision data recorded by LHCb during 
2011 at a centre-of-mass energy y/s = 7TeV. This sample corresponds to an integrated 


^The measurement from CLEO 40 is excluded from this average since S-wave contributions were not 


subtracted in the determination of the branching fraction. 
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luminosity of 1.0 fb“^. 

A test of the SM is performed by measuring eight CP-violating quantities which are 
predicted to be small in the SM. All of these are found to be compatible with the SM 
expectation, within 2a uncertainty. In addition, assuming no CP violation, the angular 
distribution of the decay products is analysed as a function of the Kn pair invariant mass 
to measure the polarisation fractions of the decay —)■ K*^K*^ as well as the magnitude 
and phase of the various S-wave amplitudes. The low polarisation of the vector-vector 
decay is conhrmed by the measurement = 0.201 ± 0.057 (stat.) ± 0.040 (syst.), and a 
large S-wave contribution is found. 

Finally, an update of the —)■ K*^K*^ branching fraction, using the 5° —)■ (f)K*^ 

decay as normalisation channel, yields —>■ = (10.8±2.1 (stat.) ±1.4 (syst.) ± 

0.6 (fd/fs)) X 10“®, in agreement with the theoretical prediction [^. This result takes into 
account the S-wave component measured for the first time through the angular analysis 
of Bg — K^tt~K~tt^ decays and supersedes the measurement in Ref. |^. 
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